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Abstract 



> 

f , High approximations of semiclassical trajectory-coherent states (TCS) and of semiclassical Green func- 

tion (in the class of semiclassically concentrated states) for the Dirac operator with anomalous Pauli 
interaction are obtained. For Schrodinger and Dirac operators trajectory-coherent representations are 
. constructed up to any precision with respect to h, h — » 0. 

00 

Introduction 

Ph. This paper is the second part of our work jl| and deals with higher approximations of semiclassical 

trajectory-coherent states (TCS) for the Dirac equation with anomalous Pauli interaction Q in an external 
electromagnetic field. The approach given in ||] is based on the concept that classical equations of motion 
can be considered as the limiting (as h — > 0) equations of motion for averaged values of corresponding 
quantum mechanical magnitudes. In physical literature this passage to the limit is justified by describing 
the dynamical states of a quantum system as wave packets localized near the position of a classical 
particle. For more detailed references, see paper JIJ and reviews j^, ^5). Such states can be obtained up 
to any power of h 1 / 2 , h — > 0, by the complex WKB-Maslov method called the complex germ theory 0]. 

In [1, 7-9], complete orthonormalized sets of asymptotic (as H — ► 0) semiclassical solutions to wave 
equations (Schrodinger, Klein-Gordon, pseudodiffcrcntial equation of the Schrodinger type for a rela- 
tivistic particle) were constructed for a scalar quantum particle in an arbitrary external field. These 
sets satisfy the coherence condition: the quantum mechanical values (averaged with respect to the func- 
tions from these sets) for the operators of coordinates and momenta are (as h — > 0) solutions of classical 
equations of motion (Newton and Lorentz respectively). 

Such solutions, called semiclassical trajectory-coherent states (TCS), generalize the well-known co- 
herent [10-14] and correlated coherent states of nonrelativistic quadratic systems [H| to the case of an 
arbitrary external field and scalar relativistic equations. 

For a relativistic quantum particle in an arbitrary external field, which is described by the Dirac 
equation, the semiclassical TCS were constructed in [16-18]. They are characterized by the fact that, 
beside the coherence condition with respect to orbital variables, they must also satisfy the following 
condition of spin coherence: in the limit as H — ► 0, the quantum mechanical averages of spin operator, 
i.e., of the Bargman polarization vector S 11 |19| , are solutions to the classical relativistic equations of spin 
motion, which are the Bargmann-Michel-Telegdi equations p(|. 
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Note that the construction of semiclassical TC-states for quantum systems is based on the semiclassi- 
cal trajectory coherent representation of the corresponding Hamiltonians. For scalar quantum mechanical 
equations, this rep resentation was implicitly used in essence in [||. For the Dirac operator, it was intro- 
duced first in |l6[ ^7|. From the viewpoint of semiclassical asymptotics, the constructed TC-states give 
the principal term of asymptotic expansion in the powers of H 1 ^ 2 , h — > 0. They satisfy the corresponding 
evolutionary quantum mechanical equations up to the functions, whose L2-norm is of the order 0(h 3 / 2 ) 
as h — > ( mod 0(h 3 / 2 )). In this paper we construct the higher approximations of the trajectory-coherent 
states of the Dirac operator in an arbitrary magnetic field, i.e., the asymptotic solutions of these equations 
up to 0{TS n+1 ^/ 2 ), where N is arbitrary integer independent of h; N = 3, 4, .. .. 

In the general situation the scheme for construction of corrections to the principal term of a semiclas- 
sical asymptotic solution is well know . Nevertheless, its realization for particular quantum equations 
with arbitrary electromagnetic potentials requires a large amount of nontrivial calculations. In particular, 
we succeeded in constructing a unitary (up to any order in H 1 / 2 ) operator in the form of an asymptotic 
series in powers of h 1 ^ 2 , H — > 0, which defines the passage to the semiclassical trajectory-coherent repre- 
sentation. 

In the semiclassical description of the behavior of a quantum particle with its spin properties taken into 
account, this allowed to construct the two-component theory, i.e., to use the space of positive-frequency 
solutions to the Dirac equations, to exclude explicitly the Hamiltonian which is a relativistic generalization 
of the Pauli operator [ pl| . The scalar part of this Hamiltonian describes quantum fluctuations of the wave 
packet near the particle position on a classical trajectory, and its "vector" part describes the interaction 
between the particle spin and the external field and the quantum fluctuations with respect to orbital 
variables. 

Note that in order to study certain specific physical problems, e.g., to consider with a given precision 
as h — > 22 , 23] the quantum effects which appear when a charge radiates spontaneously in an arbitrary 
external field, we must take into account the higher (as h — > 0) corrections to the principal term of TCS. 
In this case, to obtain the first (up to 0(h 2 )) quantum corrections to the power of radiation (which 
are uniform with respect to relativism, i.e., they can be used in any region of the particle energy), it is 
necessary to take the semiclassical TC-states satisfying the Dirac equation up to 0{h 5 / 2 ). 

The influence of corrections to the principal term of asymptotics is essential when nonadiabatic phases 
in quantum mechanics are considered pif , as well as when equations of motion of quantum averages are 
derived [1, 25-27]. 

This paper is organized as follows. Section 1 contains some relevant facts about the semiclassical 
trajectory-coherent representation (TC-representation) for the Schrodinger equation. In Section 2 we 
give the definition and investigate the simple properties of the semiclassically concentrated states of the 
Dirac equation. Higher approximations of semiclassical trajectory-coherent states \&t (a:, i, K) and of 
semiclassical Green function (x, y, t, s) (in the class of semiclassically concentrated states) for the 
Dirac equation are constructed in Sections 3 and 4. In Section 5 we show that in the class of positive- 
frequency (negative- frequency) Vh, h — ► 0, semiclassically concentrated states it can be gone to the one 
particle two-component theory whose Hamiltonian 7i^\t) is a self-ajoint operator for any order of h — > 0. 
Portion of necessary material is taken into appendices. Our paper is essentially based on the results and 
notation of Part 1 [0, some of them are not defined in this text, but only referred to as, e.g., (1.2.3), 
which means formula (2.3) in the first part Q. 



1 Semiclassical trajectory-coherent representation 
for the Schrodinger equation 

As is known, physical results in quantum theory are independent of the choice of the representation for its 
principal dynamical variables. A good choice of one or another representation often allows to simplify the 
problem or to solve it completely. To solve the problem of passing in the limit from quantum to classical 
mechanics, a new semiclassical trajectory-coherent (TC) representation was constructed in jll| up to 
mod (^(fi 1 / 2 ), h — > 0. In this representation, in the limit as h — » 0, the quantum averages for an arbitrary 
observable A%{t) with classical analog An(t) — A(p,x,t) go over into the classical observable A(p,x,t) 
translated along the trajectories of the corresponding Hamilton system (1. 1.1). In concrete problems, 
from the viewpoint of applications, it is important to know how the quantum-mechanical averages tend 
to their classical limit, i.e., the corresponding asymptotic expansions must be obtained up to any precision 
with respect to K, h — > 0. To solve this problem, we shall construct a semiclassical trajectory-coherent 
representation for the Schrodinger equation (1.0.1) up to 0{h^ N+1 ^ 2 ), N = 1,2, In Section 4, in 
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order to construct a semiclassical TC-representation of the Dirac operator in an arbitrary external field, 
we shall use the construction of a unitary operator (with a given precision in h — > 0), which defines the 
passage to such representations. 

Let us define the Hilbert space of functions depending on x and h with the following scalar product 
of two functions fi(x, t, K) and ^pi{x, t, K): 



(Pi\<P2)l% = J P*{x,t,K)ip2{x,t,h)p z ^{x,t)d x, (1.1) 

where the normalized measure density (x, t) is equal to 

p z h {$, t) = N 2 \J(t, zq)!" 1 exp Im S(x, t)}. (1.2) 

Here the complex action S(x,t) and the normalizing factor Nn are defined in (1.1.11). It is natural to 
consider the space L l n as the space of states of quantum system (1.0.2) which, as h — > 0, are localized in 
a neighborhood of the position z(t,z ) of a classical particle on the phase trajectory, since p^'(x,t) — > 
S(x — x(t)) as h — > and, in the p-representation, we have F^ p p^°(x,t) — ► S(p — p(t)), where is 
the ^-Fourier transform f2S§. The measure depends on a small parameter which implies that one can 
approximate (up to a given precision h — > 0) the smooth functions tp(x, t) by partial sums of the Taylor 
series in powers of Ax: 

N 1 

p(x,t) = ^2-d k <p{t)+R N (x,t), 

k=0 

where d k cp(t) is the fc-th term of the expansion of (f(x(t,zo) + Ax, t) into a Taylor series in powers of 
Ax, and for Risr(x,t) we have the estimate: \\Rn(x, t)\\i,t = Oih^ N+1 ^ 2 ). This estimate follows from 
the relation J (x,t)(a, Ax) N+1 dx = 0(h( N+1 ^ 2 ) which holds for any constant vector a (comp. with 
(1.4.8)). On the other hand, the space L l h contains, as its elements, the functions ip(x,t, H) which depend 
on h — ► in a singular way and are of the order O(l) as h — > in the norm of L\. For example, 



„, /(a(t),Ax)\ k 



Denote by Vn the set of polinomials in powers of h 1 / 2 Ax with coefficientss depending on t. The set 
Vh is everywhere dense in L\ and is natural domain of definition for differential operators acting in L l h . 
By 0(h a ) we shall denote the operator F(K) : L l h — > L l h for which on the set Tn the estimate holds: 
\\F(h)ip\\ L ^ — 0(h a ) as h — > uniformly in t = [0,T]. Note that in this sense the operators 

VhV, (Ax,V), ± = IL + (2(t,z ),V) (1.3) 
at ot 

are the order O(l) as H — ► 0. 

Let us define the operator t ( s N \t,H) : L\ -> L 2 (RJ), which defines (up to mod6(^ Ar+1 )/ 2 )) the 
passage to the semiclassical TC-representation by the formula: 



71=0 

"I 



4 0) (t,ft) = ^ exp ls(a?,i) , (1.4) 

oo * 

£x{t)ip{t) = M / <Mt,HWi(t)Mt)>, 

kl=o o 

where the functions S(x,t), J(x,t) and the constant iV^ are defined in (1.1.11)— (1.1.15) , and the 

operator Hi is defined in (1.4.7). 

The operator K,^\t, K) maps the space Lt into L 2 unitarily up to 0(H^ N+1 ^ 2 ) which means that 

(t { s N \t,h) Vl \tf\t,h)^ 2 ) L2 = (^Ml! L +0(ft (7V+1)/2 ). (1.5) 



3 



The Schrodinger equation (1.0.1) in the semiclassical TC-representation given the operator K.g(t, h) 



(1.4) has the form 



[tf\t, h)} ~ L { - ihd t + H}tf\t, % = n f + 0(^+ 3 )/ 2 ), 



tto = i-^){j t + [Wxp(«) + Q(*)W W (*)]V) - -^<V,W OT (t)V)}, 



(1.6) 
(1.7) 



and Q(i) is defined in (1.1.9). Thus with precision up to functions of the order 0(h^ N+3 ^ 2 ), N = 0, 1, 2, . . ., 
in the norm of the space L^, the Schrodinger equation in the semiclassical TC-representation is equivalent 
to the equation 

n <P = 0, <p&Ll (1.8) 

Equation (1.8) can be easily integrated if we note that the operator ttq admits the complete set of 
dynamical symmetries. It is easy to verify that the operators 

1 



A1 



commute with the operator ttq and satisfy the Bose commutation relations [|J : 

[A fe ,A+] =S kj , [Afe, Aj-] = [A+ At] =0. 



(1.9) 



Hence the functions H u = \H V ) = Jj {vj\) 1 / 2 (A+)' yj l form in L\ the complete orthonormalized set of 



i=i 



solutions of equation (1.8). Applying the operator K.g(t, K) to the functions \H V ), we obtain the higher 
approximations for semiclassical TCS (1.4.14) of the Schrodinger equation. 
If follows from (1.6), (1.9) thaat the operators 



[A+] {N) = (t, h)A+ (i, h)]-\ j=-, 



are an approximate symmetry operators for Schrodinger equation 



'{-thdt + n},^ 

and satisfy the following commutation relations 

~[A+r< tor 



{-ihd t + H},[A+] 



(AT) 



0(^+3)72) 



Af\ [At 



i W 



= 0(ft( Ar + 3 )/ 2 ), 

<5 Jfc + 6(ft( Ar + 3 )/ 2 ). 



For N — 2we calculate the operator K.g^\t, H) explicitly. By (1.4) and the nonrelativistic Hamiltonian 
function 



we have 
where 



H = -^-P 2 + e$, P = p--A, 
2m c 



jt {2 > {t, h) = K. { g> (t, h) [1 - i VhTd - ihn 2 - M\] , 



(1.10) 



1 

3! 
1 



(.7 + 2) 



1 r i 

-— J2 W / dT(H„\-D>+ 2 H(T)\v(T)y : j = 1,2, 
((P[,d 2 A) + (d 2 A,P[)) + |(d 3 $ - 0,d 3 A)), 



Amc 



-(d 4 $ - (p,d 4 A)) + -^—((FA,<?A) 



4! 



e 1 
2mc 3! 



(1.11) 
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f\ = -ihV + Q(t)Ax- -dA(t), /3=-$(t, z ), 

c c 

and d k A(t) denotes the function d k A{t) = ((Ax, d/dy)) A(y, i)|f=£(t, Zo ). Note that since the operators 

V'i and Ax are slef-adjoint in L l h formulas (1.11) define the self-adjoint in L l h operators D 3 T-L(t), j = 3,4. 
In particular, this implies that the operator K.g*\t, h) is unitary for N — 2 (up to 0(h 3 ^ 2 )): 

((tf(t,h)H„,)\(£f(t,h)H v )) L2 = {H v ,\(tf(t,n)) + tf(t,n)\Hu) Ll = 

= (H V '\[l - iVh(7ti - n+) - ih(n 2 - tt+) + ^tt+tti - %\ - (tt 2 )+)] \H v ) + 0(h 3/2 ) = 

= S v y + 0(h 3 ' 2 ). (1.12) 

By using (1.10) for any arbitrary operator Afi(t) — A(p, x, t) with Weyl symbol A(p, x, t) possessing a 
classical analog, we can find its explicit form (up to 0(h 3 / 2 )) in the TC-representation 

A' n = (4 2) (t, h^An^tf (t, K) + 6(h 3 ' 2 ) = 

= A(t) + t) l A(t) + ^D 2 A(t) - i^/h^D 1 A(t))n 1 - n+b l A(t)] + 6(ft 3/2 ), (1.13) 

where by & A(t), similarly (1.4.8), we denote the operator of the form 

* a ® - h + ^ v.ir.t)!^,,^,, (L14) 

Aff = -iKV + Q(t)Ax, 

and the function A(t) is the classical observable A(p, x, t), corresponding to the operator An(t), calculated 
in the point z(t, Zq) on a phase trajectory, A(t) = A(p(t, zq), x(t, zq), t). 



2 Theorem about semiclassicaly concentrated states of the Dirac 
equation 

Let us consider a relativistic particle described by the Dirac equation with anomalous Pauli interaction 

L D t> = 0, L D = -ihd t + H D , (2.1) 
where Hamiltonan Tt p has the form [B9| : 

H D = Ho + (-ih)Hi, U Q = c(a,V) + p 3 m Q c? + e$(x,t), 
Hi = te " (g ~ 2) [p 3 (S, H(x, t)) + 92 (% E(x, *))] , 



a = pi£, 1 , e = -e, 

/ o i \ / o -a \ (i o 
* = ^ i o J ' P2 = a o J ' P3 = { o -i 

Here <r = (<7i, 02, 03) are the Pauli matrices, and I are the zero and iunit 2x2 matrices, H(x,t) = 
rot A(x, t) is a magnetic field, E(x, t) — — V$(£', t) — c~ 1 d t A(x, t) is the electric field, g is the gyromagnetic 
ratio. 

Precisely as in the nonrelativistic case, we assume that the electromagnetic potentials $(:?, t) and 
A(x,t) are smooth functions of x G M 3 , i £ I 1 and as \x\ — > 00 increase uniformly in t e R 1 together 
with their derivatives not faster than any power of \x\. 

By analogy with the Schrodinger equation, we define the semiclassical concentrated for the Pauli 
equation: 

Definition 1 The state ^ will be called semiclassically concentrated on the phase trajectory z(t) = 
(p(t),x(t)) of the class CS D (N) =CS D (N,z(t),H) (* € CS D (N, z(t), h)) if 

(i) L D V = 0; 
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(ii) for the wave function ^{x,t,K) in the x-representation and for the wave function \I/(p, i, ft) in the 
p-representation there are generalized limits: 

lim |*(£,i,ft)| 2 =8{x-x(t)), (2.3) 
h— >o 

lim |* (p,t, h)\ 2 =S(p-p(t))l (2-4) 
ft— >o 

(iii) for h G [0, 1[ there exist quantum momenta ^ A°^\ \a\ + \[3\ — k, < k < N, a, (3 G Z^_, 



fc=i 

Here 

WAS) = (*|r j A c a 1 5)|*) I3 , (2.5) 
where A^ 1( ^ is defined in (1.0.13), ( | )d is the scalar product 

(*i|*2>zj = J d 3 x^+(x 1 t,h)^ 2 (x,t,H), (2.6) 

and Tj is a certain basis in the space of Dirac matrices. 
By analogy with the nonrelativistic case we have 

Theorem 1 Ifty(t) is a semiclassically concentrated state of the class CSu{z(t), N), then x(t) and p{t) 
are solutions of the classical Hamilton system with Hamiltonian \( + \p,x,t) or \(~\p,x,t), where 

A (±) (p, x, t) = e$(x, t) ± e(p, x, t), e(p, x, t) = \J ' c 2 V 2 + m%c 4 . 

Proof. 

1. We consider the spectral properties of the principal symbol of the Hamiltonian Hd, i- c -, of the 
matrix H (p, x, t) (2.2), p G R£, x G R|. The equation 

det\\H (p,x,t)- \I 4x4 \\ =0 

has [28] (see Property 2, Appendix A) two eigenvalues A of multiplicity 2 for all p G Mp, x G K^, t G M 1 : 

A^ (p, f , i) = e$(aT, t) ± £(p, f , i) 
e(p,f,t) = ^c 2 T 2 +m 2 c 4 , P = p-^A(x,t) 

We combine the eigenvectors /^(p, a;, t), j = 1, 2, corresponding to A^ (p, af, i), into 4x2 Pauli matrices 
n±(p,f,t): 

H (p,x,t)U±(p,x,t) = \ {±) (p,x,t)I\±(p,x,t), 
n + (p, f, t) = (2s(e + moc 2 ))- 1 / 2 ( ) , (2.8) 



n_ (p, f , i) = (2e(e + m c J )) 



2^-1/2/ c(a,P) 



-TOqC 2 — £ 



2. We expand the function *S? with respect to eigenvectors (on the trajectory p(t), x(t)) of the principal 
symbol of the Hamiltonian: 

* (f , t, ft) - A R (n+(t) (a?, t, h) + n_ (t) j(-) (a?, t, ft)) , (2.9) 

where n±(i) = n± (p(t), a?(i), i) , AT^ = const. The spinors are represented as follows 

T^(r t ^(aUft) exp^S^^.t.ft)) \ , . 

- I Q?\g,t,h) e Xp (iZ?\x,t,h)) J' (2 ' 10) 

where sj^, Qj^, fc = 1,2, are real functions. The condition (2.3), (2.4) implies that, without loss of 
generality, one can set 

, 3 v -1/2 

Q{ ± \x 1 t,h)=(J[a l (h)) P ( k ±} (Ct,h), Z j =Ax j /a j (h), (2.11) 
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where cij(h) are certain nonnegative functions of ft, such that 

£s^ = &^o =0 - (2 - 12) 

Since there exist momenta ^ Aq'q '' , the following limits are exist: 

lim (E { k ±] (£ t, ft) - S ; (±) (£ t, ft)) , lim (sf 5 (£ t, ft) - (£ t, ft)) , 

k, I = 1, 2. Hence, 

4 ±) (e,i,ft) = <f>o(C,<,ft) + *i ±) (e,i,ft), fc = l,2, (2.13) 

where ^^(C*, t, ft) regularly depend on ft. 

Similarly, we get for the wave function in the p-representation: 

^ = Api/bjW, lim 6,(ft) = lim^ = 0, j = M, 

and 

sf ) (^t,ft) = $o(»f,t,fi) + *i ±) (»f,t,ft), fc = l,2, (2.15) 

where $£^(77, i, ft) regularly depend on ft. 

Precisely as in the nonrelativistic case, the conditions 

\\Axj*\\ = \\Axj*\\ ~ aj (h), 
||A^|| = ||Ap,*||~&,(ft), 3 ' 



implies 



as well as 



and 



< l im , ?A m < °°' i = 1 > 2 > 3 > ( 2 - 16 ) 

ft->0 aj(h)bj(n) 



$o(£ t, ft) = *o(t, ft) + (#(*), Af) + S(£ t, ft)) / ft, 
$o(»f, *, ^ - *o(t, ft) + «#(*), Ap) + S(t?, t, ft)) / ft, 



(2.17) 



S(£,t,h)~S(if,t,h)~C(h), 

C(ft) = min{ai(ft)6i(ft), a 2 (ft)6 2 (ft), a 3 (ft)fo 3 (ft)}. 

3. We find the limit (as ft — > 0) of the averaged value of the operator with respect to the state 
(2.9), (2.10), (2.13), (2.17). By (2.1), (2.2), we get 



0^ lim ft "7 ' - (p{t),m) + A (+) (t) Hm\\j(+\x,t,h)\\ 2 + 
h^o at I h-^o 



rfgoM) 

(^limft^l^ - W),m +A ( " ) W) lim ||j(-)(x,t,ft)|| 2 . (2.18) 



The expressions in parentheses in (2.18) cannot vanish simultaneously, thus, 

$o(t,ft) - ^ (i,ft), 
t 



S (t,0) - y (^.^-^'(t))*, (2.19) 



lim || Jfr^t, ft) || = 0. 
ft— »o 



4. For definiteness, we put 

t 

S„(t,0)= f((p(t)^(t))-^ +) (t))dt, lim||j(-)(f,i,ft)|| =0. (2.20) 
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By (2.3), lim || j( + '(x, t,h)\\ = 1. Then, averaging the operator Ld with respect to the state (2.19) in the 
higher orders of H, we obtain 







(^|^+A< + )(*)-<*(*)^*)>)s:!^ (o) (t 



+ e j:j< (i) (^K(ft) 



*;=i 

(+)„cl(l)/ 



,&) + 



+ 



+ ^>a^(t,h)b k (h) 



e(i) 



+ . 



(2.21) 



Here 



{ ^°Ht,h) = \\j(+\x,t,h)\\, 
i 



!K (1) m) 



( jW (f, t, fi)) + Ax fe (x, t, h)d 3 x, 



(2.22) 



(+) 
(+) 



< (1) M) = Stfj£}J (J {+) (x,t,h)) + Ap k J^(x,t,h)d 3 x. 



We took into account that, under the conditions (2.20), the order (with respect to K) of the functions 
( ~jcr cl W, [+ ) ) cr cl(1) , and [lj<7 cl(1) is less than that of the function (2.22). Since the functions (+j<^ (1) (i, K) 



and [^jor)* *\t, K) are independent, Eq. (2.21) yields 



p=-\ { +\p,x,t), x = X^'(p,x,t). 



(2.23) 



If in (2.20) we replace (+) -> (-) and (-) -> (+), then in (2.23) we get A<+) -> A< _ ). Thus Theorem 1 
is proved. 

Theorem 2 If ^(t) is a semiclassically concerntrated state of the class <CSD(z(i),N), then the mean 
values of the quantum-mechanical spin operator, i.e., of the Bargmann polarization pseudovector S 1 * [19], 
are (in the limit as h —* 0) solutions of the classical relativistic equations of spin motion (the Bargmann- 
Mishel-Telegdi) [20]. 

Proof. For defmiteness, we assume that x(t) and p(t) satisfy the Hamilton system (2.23) and assumption 
(2.20) holds. Denote 

((t) = lim f (j(+\x,t,h)) + aj( + \x,t,h)dx, 
ft— »o J 



o"(t) = lim <*(£,*, ft) ft))!?, 
ft— >o 



(2.24) 
(2.25) 



where 



m c 



S = p 3 S + PiP. 

moc 



Then, by (2.20), we have 



7/3 



ao(t) =7<C(*),/J>, = C(t) + t^t(CW,/3) 

1+7 

The time evolution of mean values a^(t) is described by the Heisenberg equation 



a (t) = (f3,a{t)}. 



— — = hm(* 
dt n^a 



dt h 



+ -[H D ,S»]}\*). 



(2.26) 



(2.27) 



(2.28) 



Commuting the operator 5* M with Ti.jj, we get (see, for example, [29] and Appendix B): 

-<£,£)- 



^r + h [7iD > Sa] 



e (g- 2) i 
2(m c 



TOqC 



m -(p2<£,£ x P) + P3 (S, HxV)) + 0(h), 

OS i - e - 

— + -[H D ,S} = {pxE + H xS)- 

_ eo(g-2) ( ^ + ^ x g _ p3 ^j } _Z_ + -% 3 (£, #)) + O(ft). 



2m c2 



m c m c 



(2.29) 
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Substituting (2.29) into (2.28) and taking into account (A.8)-(A.14), we obtain (see Appendix B) 

§ = ir^-{E0, a)+axH) + (g ~ 2)n ' 0(0, a x H) + (a, /?} 0, E) - (a, E)). (2.30) 
at Zmocj ZrriQC 

Passing from (2.30) to the equation for a^ 1 , and taking into account (2.27), we get 

d -f = ^Ffar + { -^x»x„F™a a , (2.31) 
dr 2m c 2m c 6 

where F^ v is the tensor of electromagnetic field, r is the proper time. Equation (2.31) is the Bargmann- 
Michel-Telegdi equation. Thus the theorem is proved. 

3 Passage to the two-component theory 

Each eigenvalue X^(p,x,t) is associated with its own Hamilton system: 

p± = -\^(j)±,x±,t), x±(0,z )=x , p±(0,zq)=pq, 

x± = \f \p±,x±,t), z = (p ,x ) e Rp X , (3.1) 

Al ±5 = e$ s (x , t) - 0± , A x (x, i)>, X ( f ^ ±^ = c/ ?± 

and the corresponding variational system 

f W± = -X x %\t)W± - X x t\t)Z ± , 
Z± = X ( ±\t)W ± +X ( ±\t)Z ± . 



(3.2) 



The initial conditions for system (3.2) are chosen similarly to the scalar case (see (1.1.5)). Here we give 
the explicit form of 3 x 3-matrices in (3.2) 

Xg1(p,x,t) = \\{e($ XkXj - ± ,A XkX] )) e ^{A Xk ,X { £A Xj ))l 
X£\p,x,t) = (Xg\p,x,t)f = -- c X^(p,x,t)\\A k x .(x,t)\\, 

X$(p,x,t) = ±j\\(6 jk -t3±(3f)\\, 

calculated at the points z (t, zo) of the phase trajectory. 

From here on we restrict ourselves to solutions of "positive frequency" (corresponding to the eigenvalue 
X^(p, x,t)). Further, in functions x+, . . . related to X^ + \ we shall omit the index + in all cases where it 
does not lead to misunderstanding. The solutions corresponding to A^ - ) (p, x, t) (of "negative frequency") 
are obtained by the substitution: X^ — ► X^~\ z + (t,z ) — > z~(t, z ), H± — > n T , and (Z + ,W + ) — > 
(Z-,W.). 

Let us quantize the classical system (3.1) by the method of complex germ, i.e., with an arbitrary (but 
fixed) trajectory of a classical particle z(t, z ) we associate the complete set of functions of the form: 

W,t) = f[^=(at(t)r\0,t), (3.3) 



where 



|0, t) = N h (J(t, zo))- 1 / 2 exp [^S(x , t)] , 

t 

S(x,t) = J mt),p(t)) -A(+)(t))dt + (p(t),Af) + i(Ax,Q(t)Ax), 
o 

No(H), J(t), a + (t) arc defined in (1.2.11), (1.2.12), respectively. It is easy to see (item 2.2 Part I) that 
the functions (3.3) satisfy the equation of Schrodinger type: 

{-ihd t + X)\v,t) =0, 

X = A<+>(t) + (?(t),Ax) - (x(t),Aj>) + \ [(Ax, X x + J(t)Ax) + (Ax, A<+>(t)A# 

+(Af>,x£ ) (t)Ax) + (A^X^(t)Af)] = A(+)(t) + <5 1 A(+)(t) + ^ 2 A (+) (t). 
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The eigenvectors (2.8) of the matrix H.o(p, x, t) form an orthonormalized basis in C 4 (see Property 11 
in Appendix A): 

ir^(p,x,t)Ii±(p,x,t) = hx2, 
U+{p,x,t)U T {p,x,t) = 0, 

u k(p,x,t)Il+(p,x,t) = I 4x4 , 

fe=± 

and the system of scalar functions (3.3) is complete in P^(M 3 , C): 

oo 

(t,v'\v,t) = S v , v ; Y2,\v,t){t,v\ = \, 

\v\=0 

and 

P*(K 3 ,C) = {/,/ = cxp^(5 (t) + (p(t),^M%t,h), <f>({,t,h) e sj, (3.4) 

So(t) = /((^t),^)>-A(+)(t))dt, 
o 

where <f>(£, t, ft) is a smooth function in t £ [0, T] regularly depends on ft, and S is a Schwartz space with 
respect (eK 3 . The solution of equation (2.1) will be sought in the form 



t, ft) = ft) = (n+(t),n_(t)) 



i/(x, t, ft) 

V(z, i, ft) 



(n+w + n_v), 



(3.5) 



where the matrices (2.8) are calculated at the point z(t,z ) and the unknown two-component spinors 
U(x, t, ft) e "P£(K 3 , C 2 ) and V(x, t, ft) e "P£(K 3 , C 2 ) must be determined. 

Remark 1 Note, that if^ € P^(M 3 ,C 4 ) f<2.5) ifte solution of the Dirac equation than & is semiclassically 
concentrated states of class SC^ (z(t), oo) . 

We substitute the function (3.5) into (2.1) and expand the obtained expressions with respect to the 
eigenvectors (2.8) taking into account the relations (see Appendix A) 



Pl (a,V)U ± (t) = ±0,V)n ± (t)+n T (t)((aJ)^^ - (a,V) 

V 1+7 

P2 (E, E)Il ± (i) = -n± (i) (a, ^x£)T zn T (t) f (a, [3) + 7" 1 <<?, £) 

V 1 + 7 



p 3 (X,H)Il±(t) = T n ± (t)^(a,0}^-^ - (a,H)j + (0,H)Tl T (t) 
Here cd — x(t, zq), 7 _1 = \J\ — 2 . As the result we get 



(-ihd t + n)^ = IL, 



+ h{a,0x8) e(g-2)h 



2 1 + 7- 1 4m c 



jft<9 t + A< + > (t) + e A$ + c(/3, AP) + 

((a, i?) - (<?, ^xl)- (a, |)M)] W(f] t) + 



e( 9 -2)ft/._ 3 . 



4moc 



n_ 



a, /3) + h - 1 E) + 0,H))]V (x, t) } 

ihd t + A(")(t) + eA$ - c0, A?) + ^f'f X f - 

2 1 + 7 
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4m c V v ' N ' 1 + 7- 

- £ %^^ ( - * ^' ^ rrr^ " ^" 1 <5? ' ^ + ( ^ # >) ] *) } ' ( 3 - 6 ) 

where A$ = &(x, t) - ${x(t, z ), t), AV = V - P(t). 

We transform the obtained equation by expanding the expressions in this equation into a Taylor series 
in operators Ax and Ap and taking into account that in the class of functions V\ we have 

Ax = 6(Vh), Af=6{Vh), (-ihd t + \) = 6(H), 
as h -> (see (1.4.4)). Denote V\ = Ap — {e/c)d l A{t), 

k=2 ' u fe=0 ' 

- i^ 1 (a,d k E) + 0,d k H)^} + 6(ft< Ar + 3 )/ 2 ), (3.7) 
^ w - 2e(t)-2c<#Pi> + + A)- l(^ 1 ,AWf 1 ) + 

- eJ ^£±{-<*J + (3-8) 

fc— 

Then equation (3.6) takes the form 

{-ihd t + = n+ | [ - thdt + a - 1(^1, a(+'A) + 

' k=3 
fc— 

+ )+ v ^ (x, t) J + n_ [jw w v w (x, t) + #">w w (x, t)] + 6(n( JV + 3 )/ 2 ) . 

Since the eigenvectors of the principal symbol of the Hamiltonian Hd, which define the matrices 11+ (i) 
and II_(i), are linearly independent in C 4 , then, by (3.5), we get 

vW(f,i) = -{^Wj-^WyW^i), (3.9) 
H^ ( + F (A,) )W w (i,i) = 0, (3.10) 



fc! 

fe=3 

JV 



2 1 4m„c f-; fc! V 7 x 1 



k=0 



2 4m c 

_ (aJ) & dkH h - {RWy+iMwyLRW + 6{^ N +^/ 2 ). (3.11) 
1 + 7 > 

The operator (.M^) -1 inverse to the operator (3.8) will be found from (1.4.10) where we put B = M — 2s, 
A = 2s. Then 

^-g (Ar) - -{jCfW}- 1 ^) = 1 V ft fc / 2 Q fe + 0(ft< Ar + 3 )/ 2 ), (3.12) 

/c=i 



n 



where 



VfrQi = (3.13) 



and for k > 2 



- i(ff,fi { ^~[f - ^(ff, d k - 2 E)) - ^0,A)Qk-i + [e0, d 2 A) + (-iW t + A) 



+ (a,/3 x %)})] — ^— Q fe _ 2 + £ [^(d»* + 03,<fM» + 



n— 3 



4m c (n-2)!V 7 x 7 x I + 7- 1 



Qk—n • 



2e 

By taking into account (3.8),. . . , (3.11), the operator can be represented in the form 

where 

Fo = A+ 2 l+ 7 -i ~2TK + (^H(t))) + 

+ e -^T^ d{t)) <*>0 xS ®) (*>fiT^), (3-14) 

and the highest terms have the form = 6{K) as ft — > 0. 
Thus we have 

fW(t) = (n + (t) + ±n_(*)QW(t)) (3.15) 

and the problem of constructing asymptotic (up to 0(h^ N+1 ^ 2 )) positive-frequency solutions ( + )ty(x, t, ft) 
to the Dirac equation (2.1) is reduced to solving the equation (3.10) with respect to the two-component 
spinor U^ N \x, t). As a result we have 

t, ft) = *(2,t,h) =f {N \t)U {N) (x,t). 

4 Semiclassical trajectory-coherent representation of the Dirac 
equation with anomalous Pauli interaction 

Now let us consider the construction of asymptotic (as ft — > 0) solutions of equations (3.10) in the 
two-component theory: 

(-ihdt +F Q + VhP[ N) )U {N \x,t, ft) = 0. (4.1) 

In the first approximation as ft — > (neglecting the operator \fhF\ = (9(ft 3 / 2 ) in (3.10)) we have the 
following equation for the spinor U^\t): 

{-ihd t + F )U {0) {x,t,K) = (4.2) 

and hence (x, t , ft) = (j)(x,t,h)u(t), 

{ l Jt + Ye[ {l + ~ gi){ °> " {t)) ( TTT T + ~ giWl ^ x ^ (<)> ~ 
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^-(a,0)]}u = O, (4.3) 



1 + 7 

(-ihd t + \)4>(x,t,h) = 0, 
where g — (g — 2)/2. The derivation of (4.3) was based on the following relation [30]: 

'-* BC -* -* -» -* -* -> 

= —{E + x H - 0{0,E)). 

£ 

Thus the spinor properties of electrons in the semiclassical limit as h — > are determined by complex 
solutions of the linear system (4.3). Assume that at the initial time the spinor u(t) satisfies the following 
condition pl| : 

(*,i)u(o,o = C«(Q.O, C = ±i- (4.4) 

By this assumption we fix that at t = the spin of a particle is directed along an arbitrary unit vector 
£ e K 3 . Then at any instant of time the solutions u(t, (,) of the Cauchy problem (4.3), (4.4) form an 
orthonormalized basis in C 2 : u+ (t, £')it(i,£) = $(.C' an d, hence, the functions 

w$(f,t,79 = u(t,OM = kC,*> 

form the complete orthonormalized set of solutions of the equation (4.2) 

{tX',v'\v,^t)=5 v , v .5 u .. 

The precision of constructed solutions can be estimated similarly to the case of the Schrodinger equation 
(see Theorem 2 in [1]). 

Let us construct the following (in h — > 0) approximations for equation (4.1). In contrast to the 
nonrelativistic case, the operator of "perturbation" ^/Hf[ N ^ is not self-adjointfj] in L2(i? 3 ,C 2 ) and hence 
the function U^ N \x, t, K) cannot be considered as the wave function of electron. However, to solve 
equation (4.1) by the methods of the theory of perturbations and then to go over to the semiclassical 
TC-representation, it is sufficient that only the operator Fq be self-adjoint. By using formula (1.4.10) 
with operators 

1 = a + ^ ' S 'Tn F '"' 

and taking into account that 

oo 1 

i"V(*)= E E I^C,*) / <*r<T,C,i/|0(T)>, 
M=0C=±i o 

we obtain the solution of equation (4.1) with precision up to 0{h^ N+1 ^ 2 ): 



N 



OO „ 

p[ N) m = E E w,c,t) dT{TX',AF[ N) \m)- 

k'|=0C'=±l n 



k'|=0C 

Now let us construct the semiclassical TC-representation for the Dirac equation (2.1) following the 
scheme of constructing the TC-representation for the Schrodinger operator. We introduce the Hilbert 
space of vector-functions 

I4(K 3 ,C 2 ) 

with scalar product 

{Vi\<P2)l\ = J d 3 xp z h °(x,t)tp+(x,t,h)ip 2 (x,t,h) 7 
where the density of measure p^°(x,t) was defined in (1.2). 

1 A similar situation takes place for the highest nonrelativistic approximations |32|. 



13 



We define the operator JC^\t, K) : L\ — > ^(R 3 , C 4 ), which defines the passage to the semiclassical 
TC-representation, by the formula: 

(t, h)<p = fWfWjt^ (t, h)ip, y £ 4(R 3 , C 2 ), (4.6) 

where the operator T( N 1 is defined in (3.15), K.p(t,h) — K,^ (t, K), and K.g(t,h) is defined in (1.4), in 
which the symbol TL{p,x,t) must be replaced by the relativistic Hamiltonian function X^ +s> (p, x, t) . 

The operator K^\t, H) (with precision up to 0(h^ N+1 ' > ^ 2 )) unitarily maps the space L^(R 3 ,C 2 ) into 
the space L2(K 3 , C 2 ), this means that 

(^\t^) Vl \^\t,h)^ 2 ) L2 = (V>1|V2> £) > +0(^ +1 )/ 2 ). 

By direct calculations it is easy to verify that in the semiclassical TC-representation the Dirac equation 
(2.1) takes the form: 

ft™ (t, h))+{~ihd t + n D )K,^ ] (t, h)<p = [tt + (<?, A(t, z )}} <p + o(^ Ar+3 )/ 2 ), (4.7) 

where the operator ttq is given by the formula (1.7), in which the symbol H(p, x, t) must be replaced by 
the relativistic Hamiltonian function A^ + ) (2.7), and the vector T><y{t, z ) is equal to 

V (t, z ) = ^ f(l + §7)H(t) - (— L_ + ~ gi )p x E(t) - —P—00, H(t))} , (4.8) 
7 L Vl+7 1 / 1 + 7 1 J 

where fj,o — heo/(2moc) is Bohr magneton. 

Thus, with precision up to (J(h^ N+3 ^ 2 ), the Dirac equation in the semiclassical representation takes 
the form: 

{n + (a,'D (t,zo)))v = 0. (4.9) 
Precisely as in (1.8), equation (4.9) admits the complete orthonormalized set of solutions of the form: 

|^C)= M («)n-^(A + )M, (4.10) 



k=l 

where the operators of creation are given by formulas (1.9), in which the complex vectors Wj(t) and 
Zj(t) are solutions of the variational system (3.2). 

We return to the initial (Schrodinger) representation of the Dirac equation (2.1) and, taking into 
account (4.6) and (4.9), obtain the complete orthonormalized set of semiclassical (up to 0(H^ N+1 ^ 2 )) 
trajectory-coherent states of electron: 

*2 ) (j?,t > ft)=^ ) (t,ft)iH 1/ ,o. (4.H) 

The negative-frequency semiclassical TC-states ( _ )\f r £, ^ are given by the same formula where we change: 
A<+> — A(->, n±(t)— »IL F (t), z+{t,zo)—>z-(t,zo), 

(w + (t),z + (t)) — > (#_(«),£_(<)). 



Remark 2 The relations (4.6)-(4.11) allow us to obtain the Green function for Dirac equation in semi- 
classical trajectory-coherent approximation. Analogously (1.6.4) for the positive-frequency part of a kernel 
of evolution operator for the equation (2.1) 

G W (x, y,t,s) = T W (t)f^ (t)G ( ff (x, y, t, s) (,)) + (f W („)) + , (4.12) 

where 

G$\x,y,t,s)=G(°Xx,y,t,s) £ u(t, C)u + (s, C), (4.13) 

C=±i 

u(t, Q is a solution of equation (4.3) with initial condition (4.4), G^°\x,y,t, s) was defined in (I.Al.l) 
for Ti.o = A. The operators T( N > and J-^ N ^ were defined in (3.15) and (4.5), respectively. Analogously 
to the scalar case, the Green function allows to obtain a solution of Cauchy problem for Dirac equation 
(2.1) only in the class of semiclassical trajectory- coherent states. Nevertheless the Green function (4.12) 
can be usefull for calculation of concrete physical effects |29, |33[]. 
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5 The relativistic analog of the Pauli equation 

The semiclassical description of a quantum particle with its spin properties taken into account allows 
(with an arbitrary precision in h — * 0) to exclude the interference between the positive-frequency and 
negative-frequency states (the Schrodinger "zitterbewegung" J2{J), i.e., in the subspace of positive- 
frequency (negative-frequency) states it allows to go over to the one-particle two-component theory, 
whose Hamiltonian is a self-adjoint operator for any order of H — > 0. 

The unitary operator TW ( mod 0(h^ N+3 ^ 2 )) of transition to the two-component theory will be sought 
in the form: 

f W = f WbW + 6(^+ 3 )/ 2 ), (5.1) 
where the operator is defined in (3.15) and the operator B( N ) is defined by the condition that the 
operator 

(+)^(iv) is S elf-adjoint: 
In this case the spinor 

can be considered as the wave function of the one-particle problem: 

(T^)+(-ihd t + # d )tWz?M = (-ihd t + (+W N) )U^ + d(^ N + 3 ^ 2 ). (5.3) 
Let us consider the construction of the operator B( N ) in the case N = 2 more precisely. For this 

" (2) 

purpose, we represent the not self-adjoint part of the operator Ff as follows: 

l^Q x {-iKdt + h)^Qi = -{(-tfifit + F )(y £ Qi) 2 + (^Qi) 2 H^ t + F )} + 

+ ^{[(^Qi),Hm + F )]_^Qi + ^Qi [(-m + F ), 

= AB-BA. 

If in Eq. (5.1) we choose the operator B^ as 

5 (2) = l-^Qi) 2 , (5.4) 

then it can be easily verified that the conditions imposed on the operator T( N 1 are satisfied. In this case 
the spinor U v £ takes the form: 

W$(5Uft) - {B^y 1 U%{x 1 t 1 h) = 

= [l - iV0i - ihT 2 - hF 2 ]\isX,t) = F {i) \v,Q,t), (5.5) 



where 



M=o,c=±i o 



M=o,C'=±i 



- ^(-ifiar + f )(1q 1 ) 2 - J(^Q0 2 (-^ + Fo)} \<p(t)). 

Thus the functions (5.5) form the complete orthonormalized set of solutions of equation (5.3): 

which can be considered as the relativistic (up to 0(H 5 ^ 2 )) generalization of the Pauli equation. In the 

case N > 2 the operator B^ (7i.( N \ T^O respectively) can be obtained in the same way as the operator 
(5.4). 
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~ (2) 

Let us explicitly calculate the operator fC D (t, ft) which defines the passage to the semiclassical 
trajectory-coherent representation. By (4.6), (5.4), we get: 



(i - ^Qi) 2 ) [i - - A - hT 2 ]\t^ ] (t, ft), 



(5.6) 



where K^\t,H) is defined in (4.6). Since the operator is unitary, the operator (5.6) is also unitary, 



i.e.: 



{{t%\tMHv<'m^\tMH„c,)))L 2 = 

= {CH v '\H v QlI + 0(fi 3/2 ) - 8 v y8 u > + 0(ft 3 / 2 ). 

Let A t : L2(R 3 ,C 4 ) — > L2(K 3 ,C 4 ) be a unitary operator. Then the corresponding operator in the 
two-component theory can be always represented in the form: 

X+ = {T^j+AtT^ + 0(ft 3 /2) = a+(a,A}+ 6(ft 3 / 2 ), 



where a and A are self-adjoint (in L 2 ) operators with symbols a[x,p,t) and A(x,p, t) respectively. We 
find the explicit form of the operator A t (ft) in the semiclassical TC-representation: 

A' = (K^(t,H)) + A t }C^(t,h) = {t^\tM {+ H^ 2) ) (+) KF {2) t%\t,h) = 
1 + D 1 + i^ 2 ) (a(f) + {a, A(t))) - iVh^D^att) + (a, A(t)))ni - 



tc+D 1 (a(t) + (a, Alt)) + K(A(t), [tt+ cttti - 

i(a7c 2 - tt+ct) - (ttt 2 - (tt+) 2 ct]) - iVh(A(t), (otti - tt+g?)) + 0(ft 3/2 ), 



(5.7) 



where 



nj = (K^\t,h)) ^-/C^M), j = l,2, 

and WA{t) are defined (1.14). 

In particular, in the semiclassical TC-representation, the operators of momentum p = — iftV, coordi- 
nates x, and spin 2 S — ha/2 have the form: 



k{t,h) = {ic i ° ) (t,h)}- 1 {T^}- 1 xT^ic^(t,h) = 

= x(t, z ) + Ax ~ iVh(AxTri - Ax) + 6{h i/2 ), 

P(t,h) = {^(^ft)}" 1 ^ 2 )}- 1 ^)^ ^ 

= p(t, zq) + Ap- iVh{ApTt x - TT+Af) + 6(ft 3 / 2 ), 



(5.8) 

(5.9) 
(5.10) 



Let us write the explicit expression for the Hamiltonian of the two-component theory 7iW (5-2), in which 
the operators of the order 0(ft 3 / 2 ) are taken into account, and the expressions for quantum averages of the 
principal observables in the theory, namely, for operators of coordinates, momenta and spin are calculated 
for one-particle semiclassical TC-states of electron |i? V! C) (4.10): 

^ = (l + ^ + ^ 2 + ^ 3 )A«(t) + 

+ ^{ [(1 + 97M H(t)) + (— + g 7 M P x E{t)) + 
I L 1 + 7 



2c 



1 + 7 1 JeL 1+t 1 J 



2 In the initial (Schrodinger) representation, the spin operator S = (h/2)cr corresponds (up to 0(?i 3 / 2 )) to the three- 
dimensional unit vector of spin (h/2)a° = (h/2){p 3 f, + p 1 (c/e)V - cp 3 "P(S, V)/[e(e + m c 2 )}}} (see |||). 
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91 



(a, (3y.dE)- (a, dH) + (a, (3) 



{f3, dH) 



i + r 



1+7 



1 + T 



>} 



= (l + 5 1 + i^ 2 + i<5 3 ) A« (t) + (a, D (t)) + (a, [D x (t)Ax + D p (t)A0\); (5.11) 



-7" 

+7-^(*)l^AD- 



1 + 7" 



e 



. 1 + 7 1 



_ c 

57- 

e 



1 + 7 1 



eLl + 7 1 J Vc / 



(- V-A(t))] 



_ c 
37- 

£ 



7 

(&£(<)) 



, r 4 



e r il 



d l A(t)),0) 



1 + 7" 1 
-l ; 



c 



1+7 



'At))) 



1 + 7" 



(-.3)^f^+7- 1 (^^) 



37 



(a, f3xdE) - (a, dH) + (a, 0) 



(P,dH) 



1+T 

= x{t, zo)5 u > - iVhiC, H u \(Ax7n - n+AS)\H v , Q + 0(ft 2 ); 

P(t,(,C) = (C,H v \P(t,h)\H v ,() = 

= Rt,z )5 u ,-iVh{(',H u \(Atfw 1 -n+Atf)\H v X)+d(ti 2 ); 

s(t,(,C) = {t',H v \§{t,h)\H v ,Q = CO- 



1 + 7- 1 



(5.12) 

(5.13) 
(5.14) 



Here = -thV + Q(t)Ax, g = \(g-2), ff(t,C,C) =U+{tX')SU{tX) is the solution of the Bargmann- 
Michel-Telegdi equation [M in the rest system 

2 

f] = -ifx V {t,z ), 

where the vector T>o(t, Zq) is defined in (4.8). 

By using the function |if„,£) (4.10), it is easy to calculate (with the same precision 0(h 2 ), h — > 0) 
the correlation matrix |Tq] , which describes quantum fluctuations of dynamical variables Xj(t, K), Pj(t, K) 
with respect to their average values (5.12), (5.14), a XuXj , u Pi , Pj , and their correlation a PijXj , i,j — 1,2,3. 
Here 

o AB = 1 -{(AB + BA))-{A){B). 

We have 



o xx = -[C(t)D- 1 C+(t) + C*(t)D- 1 C t (t)]. j 
o pp = J [B{t)D- l B + {t) + B* {t)D- l B\t)] , 
a px = \[B{t)D- x C+{t)+B\t)n- y C\t% 



(5.15) 



where C(t) and B(t) are defined in (1.1.3), and D v 1 denotes the matrix 

2u k + 1 «. 



D- 1 



Im b k 
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Differentiating the relations (5.12)— (5.15) with respect to t and writing the right-hand sides of obtained 
relations in terms of variables X, P, if, A 2 , we get (up to 0(ft 3/<2 )) 



z = Jd z \^ (z, t) + ^(d z ,A 2 d z )Jd z \ (+) (z, t) + JD z (z, t)fj 
A 2 = J(\^Y> z (z,t)A 2 - A 2 (A«)' 2 ' 2 (z,i)J, 



A 2 = A 2 , 



(5.16) 



where 



z = (P(t,h),x(t,h)), 



A, 



' pp u px 



D z rj = (D p ff,D x ff), f3=~\f(z,t), 

vector D is defined in (4.8), and 3 x 3 matrices D p , D x - in (5.11), P(t,H), X(t,h) in (5.12), (5.13). 
The system (5.16) is a closed system of ordinary differential equations for quantum averages in Dirac 
theory. The problem of correspondence between the obtained equations and the well-known equations for 
a particle with spin (e.g., Frenkel etc. |36, 3T\ ) requires a more detail consideration. The initial conditions 
for system (5.16) is choosen as follows: 



lt=o 



zo, 



Lt=o 



A 



lt=0 9 ^ 9 I 

'l-(£,ky 



(5.17) 



where k — (0, 0, 1), £ and Aq are defined in (4.4) and (1.5.23), respectively. Let us write the system (5. If 
in the form 



z = Jd z \ {+) (z,t) + ^(d z ,AA 2 A + d z )Jd z \ i+ ' ) (z,t) + JD z (z,t)U + aU, 
A = J(A^ 



d 1 _ - 
l dt + H {<T > V ° 



zz (z, t)A, 
w = 0. 



(5.18) 



The initial conditions for spinor U are defined in (4.4) and 

A I 

where 3x3 matrices Bq and Co satisfy the conditions (1.1.4) and (1.1.5). 





' Bo 


B* \ 


\t=0 


v Co 


Co J 



6 Conclusions 

Let us briefly review our consideration. To our opinion, from the physical point of view, the "passage 
to the classics" in quantum mechanics must inevitably cause the introduction of the notion of classical 
trajectory, which primarily (on the postulational level) is alien to quantum mechanics and must be taken 
from the outside. It is essential that one can construct the complete set of approximate solutions of 
the Dirac equation with the property: as h — > 0, the quantum-mechanical average values of coordinates 
and momenta are general solutions of classical Hamiltonian equations. In literature |38| it was repeated 
over and over that this possibility is not obvious. Usually (see, for example, |54|, p. 69-70) the authors 
restrict themselves to a verbal formulation of conditions imposed on the relativistic wave function of 
semiclassical type, having implicitly in mind that these conditions can be met without any essential 
difficulties. However, as a rule, such states were not presented. Our consideration shows that in order to 
construct such states explicitly, one can use the method of complex germ Q . 

Not only the principal possibility of obtaining approximate (in H — > 0) solutions of the Dirac equations 
(up to 0{h N / 2 )) for any N is shown, but also a constructive method for obtaining the corresponding higher 
approximations is given. It is essential that the expansion into an asymptotic series with respect to h 
contains half- integer powers of ft, i.e., there are series in y/h (in contrast to the standard semiclassical 
expansion in ft, ft — > 0, given in all manuals of quantum mechanics). 
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The possibility to construct such states (called trajectory-coherent states) explicitly leads to nontrivial 
conclusions. For example, one can obtain, in the most transparent and natural way, the "classical 
equations of motion" for average values of quantities, which cannot be exactly well-defined in the classical 
sense (e.g., spin). It seems natural that, for the classical vector of spin, we obtain the Bargman-Michel- 
Telegdi equation. However, in the case of arbitrary (and not only homogeneous) electromagnetic fields, 
one can explicitly show, which is nontrivial, that in this equations the fields must be taken on classical 
trajectories. Up to now (see, for example, P4] ), this fact was justified only by verbal arguments. 

In conclusion, we note that, by using the Maslov complex canonical operator, one can go over to the 
approximate trajectory-coherent representation, in which the classical trajectory is considered already 
for the approximate Hamiltonian of quantum theory. 
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Appendix A. Properties of matrices H±(t) and \^\t) 

Property 1 The following relations hold: 



[A<i»w] 1/2 



W) 



1+7 



_i u « 



3x3 



[A«(t)]- 1/2 = - 



W) 



1 +7~ 



3x3 



fc||3x3, 



(A.l) 

(A.2) 
(A.3) 



whe 



Proof. We can verify directly that 

3 



E 



pith 



e(t) \ 1 + 7 



(5 jk - Wk) = A« fc (t), f3 2 = l-r 



E 



Mi 



- Si 



Sik + 



7/3 2 



1 + 7" 1 1 1+7" 



T - 7 + 1 - fa 



1 + 7 -1 



3 



c 2 /: Q Q s e(t) 

z 

as was to be proved 



iii) E ^y(*« - PiPi)-y-(8«< + 7 2 A/9fe) - Sjk + PjMj 2 - 1 - 7 2 /? 2 ) - S jk 
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Property 2 The following relation holds: 

H (t)TL±(t) = X (± Ht)Tl ± (t), 
where Ho(t), X^^t), and n ± (t) are defined in (2.2), (2.7), and (2.8), respectively. 
Proof. Since (3 — cV/e, we get 



(A.4) 



n (t)n + (t) = 



e + mtyc 



e$ + m c 2 c{a,V) 
Sy/2 + 2-f- 1 V c(ff,P) e$-m c 2 

1 / (e$ + m c 2 )(e + m c 2 ) + c 2 (ff,P) 2 

e^/2 + 2'y- 1 \ c(ff,P)(e + m c 2 + e$ - m c 2 ) 



(e + e$)n + (i). 



Here we took into account that c 2 (ff,P) 2 = (PV 2 = e 2 — m\P . For the lower index, relation (A.4) can 
be proved similarly. 



Property 3 The following relation holds: 



0, V) 



Pl (X,V)Il ± (t) = (a,V)Tl ± (t) = ±(/3,V)n±(t) + U T (t) I (ct,/?)^—^- - (ff,V) 



= ±0,r)u ± (t) + ^-n T (t)(a, (\£Ht)) 1/2 v). 



(A.5) 



Proof. Actually, 

(a,T)Tl + (t) 



y/2 + 27- 1 V (ff,V)(l + 1 ~ 1 ) ) ' 
(ff,V)(ffJ) = 0,V)+i(ff,Vx(3) = [1 + (1 + 7- 1 ) - (1 + -y" 1 )]^, V) - 
-i(ff,Px p) = (1 + 1 - 1 )0,V) + [(1-^ 1 )0 1 P) - (ffJ)(ff,V)] 
[3 2 1 



(1-7- 1 ) 



-p,PY 



1+7 -1 1+7" 
= (l+^)[0,P)] + (aJ) 

(l+ 1 - 1 )(a 1 V) = (l+ 1 - 1 )(cr 1 V)- (IP) iff. 1) + { lP)(ff. 7) 
= (ffJ)[0,P)}-(l+ 1 - 1 ) 



(ffJ)^-(ff,P) 



P- OrT^ 



1 + 7" 

Thus, we obtain the first relation in (A.5) for upper indices. Similarly 

1 



(S,V)U. 



-(a,P)(l+ 7 - 1 ) 
P,P)pJ) 



V / 2(l+7- 1 ) \ 

-(a,P)(l+ 7 - 1 ) = 0,P)(P,ff) + 0,P)(ffJ) -(l + l- 1 )(ff,P) = 



(ffJ)^-(ff,P) 



= (a,/3)[-(/3,7')] + (l+7- 1 ) 

/3> = (1 + 7" 1 )(/3,^> + (1~" 7 * ^) - ^ 7 ) P- V) 

= (-l-l- 1 )[-0,P)] + {ffJ) 



1 + 7" 



as was to be proved. The second relation in (A.5) immediately follows from (A.l). 



Remark 3 If in the matrices U± (t) = H± (p, x, t) 
in the obtained expressions we must put j3 = cP/e. 



trajectory x(t), p(t) is not classical, then 



p—p(t),x—x(t) 
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Property 4 The following relation holds: 

n ./) = n (/) 



2 I+7- 1 



1+7" 



n±(t){?,(AW(t))- 1/2 4 



Proof. The latter relation in Eq. (A. 6) is a direct consequence of Eq. (A. 2). We consider 



dt 
d 1 



1 \ 

(<?,/?> 



dt 



1 dr^Jdt _ 1 1 (fry- 1 ( p 2 _i 

2 yT+T 1 ^ 2yr^^r^r I i+ 7 -i +7 



d (<?,/?) 



Similarly, 



ii) 



1 



2^1 + 7- 1 
dj- 1 



jP dry- 1 ] ^jd-y- 1 

I + 7- 1 dt 



2>^> 



1 



7(A/3) 



\A + 7 



= (1+7- 1 ) 



dt 

i (a, (3 x /3) 
2(l + 7" 1 ) 



2V1 + 7- 1 
1 

2^1 + 7- 1 



1 + 7 1 



u 1 - do-- 1 

2(a,/3)- T - rT 3 T (a,/3)^ r 



(<?,/?> 



1 + 7 

/3 2 



1 + r 



1 + 7" 1 



= \V 2 (aj) - (aj)(pj) = -(a J x x /?)) = x $| = 



1 



2^1 + 7- 
1 

7l + 7^ 



— — ^ + (1 + 7~ 1 )(^ 



- *\ d 7 _1 



2 I + 7- 1 



+ 2(1 + 7^) 



( ff '^ 1+7 -l +(^P) 



d 1 

dt w V2 



( ± (aj) \ 
dt 

_d_ 
V ~di 



d (<?,/?) 



dt y/l + j- 



d 
~dt 



1 



2 V /IT7 _T 

+ (l + 7" 1 ) 
1 



' 1 + 7" 1 



+ 



2 V /IT7 _T 



1 +7" 



-(1 + 7^) 



1 1 + 7- 1 



+ 



as was to be proved. 



21 



Property 5 The following relation holds: 

p s (t,H)u ± (t) = T n±(t) 



g x (0,H) ~ 



+ n T (t)(/3 7J ff) 



(A.7) 



Proof. The latter relation in (A.7) follows directly from Eq. (A.l). Precisely as in Property 3, we get 



p 3 (^H)U + (t) = 



(a,i?)(l+ 7 - 1 ) = -(l+ 7 - 1 ) {aJ)-^^-{5,H) + (a, 0)(0, H); 
-(a,H)(a,0) = -(a,0) (a, 0) ^3 - (a, H ) - (1 + 7 - 1 )(/3, 



ii) 



p 3 <s,H)n-(t) 



1 ( (^,H)(a,0) 



Further, the proof coincides with that of Property 3. 
Property 6 The following relation holds: 

p 2 {E,E)U ± (t) = Il ± (t)(-(ajx E)) T m T (t) 



1+7 1 



-U ± (t)(aJxE)± ll - 1 ^n T (t)(<j,(Xi+\t))-^E). 



/e{t) 



Proof. The latter equality in (A. 8) follows from (A. 2). We consider the relations 



p 2 (X,E)Il + (t) = 



i ( -(Z,E)(a,0) 



(A.8) 



ii) 



P2 (£,£>n-(t) 



-i(a,E)(a,0) = -i(0,E) + (a,E x 0) = 

= + 7 - 1 )(A £) + * 7 ~ 1 £> + i(ff, 0xE))+(l + 7 - 1 )<ff, £x(3> 



(l + T- 1 )!-^^^]-^,^) 



1+7" 



-+ 7 - 1 (a,i?) 



i<<?, + 7 - X ) = z fl±Xi + 0n (<?, E) 



-i((ff, 0)(a, E) - (ff, E)0 2 ) + i(a, 0) (<?, E) + ±±1-L(ff, E) 



= (aJ)[-(aJxE)]+i{l + 1 - 1 ) 
as was to be proved. 
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Property 7 The following relation holds: 



(E,s>n±(i) - u ± (t) 



7- 1 <-,S) + -^(a,/3) 
1 + 7 



±n T (t)[i(aJx S)} = 



= -^U ± (t){a,(X^(t))-^S)±iU T (t){aJxS). (A.9) 



Proof. Since 

1 ( (<?,£>(! + 7" 1 ) 



(E,5)n + (t) = 



the further proof coincides with that of Property 6. 
Property 8 The following relation holds: 

p 3 n ± (t)=± 1 - 1 Il ± (t)+Il T (t)(aJ). (A.10) 

Proof. 

i} p3U+it) = vm^) ( 

(1 + 7" 1 ) = 7" 1 (1 + 7" 1 ) + /? 2 = (1 + 7~ 1 )7~ 1 + «<?, /3)) 2 ; 
-(a,/?) = 7 - 1 (a,/3)-(l+ 7 - 1 )(a,/3), 

as was to be proved. 

ii) p 3 n-(*)= , 1 ( 

^2(1+7-1) V i + 
Proof is similar to the preceding one. 
Property 9 The following relation holds: 

P2 Tl±(t) = ^Tl T (t). (A.ll) 
Proof directly follows from the definition of matrices II±(i). 
Property 10 The following relation holds: 

Pl n ± (t) = ±n ± (t)(aJ)- 1 - 1 n T (t). (A.12) 

Proof. Multiplying the left- and right-hand sides of (A. 10) by ip 2 and transforming the right-hand side 
of the obtained expression according to (A.ll), we get (A.12). 

Property 11 For matrices H±(t), the orthonormality and completeness relations hold: 

I4(*)n±(t) = I 2X 2, I4(*)IM*) - 2x2 ; (A.13) 

n+(i)n* + (t) + n_(i)n*_(t) = i 4x4 . (A.u) 

Proof. By the straightforward verification, we get 

n 4 ± Wn±(t) = 2(1+ 1 7 _ 1) [(i + 7- 1 ) 2 + ((^/?)) 2 ] = 

! -(l + 2 7 - 1 + 7 - 2 +/3 2 )I 2x2 =I 2X2 ; 



2(1 + 7 -i) 



l4(*)IMi) = 9n f 1 ZT7 [<^/3)(l +7" 1 ) - (1 + 7- 1 )(g,/3)] = 0; 
2(1+7 ) 

nmn'm + nmn'm 1 f f (1 + 7- 1 ) 2 ( ( ?,/3)(i + 7 - 1 ) V 

n + (i)n + (i) + n_(i)n_(t) - . 2 1 + 



+ -(^xi+7- 1 ) d+7- 1 ) 2 J r 4x41 



as was to be proved. 
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Appendix B. The Heisenberg equations for the polarization op- 
erator 

Let us write the Heisenberg equation for the polarization operator: 

±(9\S^) D = (*\^\*)n + l j_(mn D , S,m D , (B.l) 

where 

S^ = (S ,'S), So = —&?), g=p 3 Z+— Pl P, (B.2) 

TUqC UIqC 



Hd =r Ho — ihHi , Hq = c(a, V) + p 3 m c 2 + e$ 



For this purpose we need following auxiliary propositions: 
Lemma 1 The following relation holds: 



Proof. Actually, 



dt m c 2 ' dt 

[Ho,§ ]- = — (S,V$>. 
m c 



Since E = — V$ — the property is proved. 
Lemma 2 The following relation holds: 



Proof. Actually 



dt moc 2 ^ 1 dt ' 

[7io,p 3 £]_ = c{{ip 2 ){P-i^xP-ip 2 {P + itxP)} = 
= -2icp 2 P, 
1 - ^ 77?P 1 - i. % 

-[Wo.PiP] = 2 ^P2^ + piV$ + — [<£,p),'P]_. 



(B.3) 



^ + |[Ho,5 ]_ = — (E,£>. (B.4) 
at ft m c 



aS'o e - OA. 

"V 2 "' "qT/- 



^ + 1[W ,<?]_ = MlfixS). (B.5) 

at n, m c 



dS e dA 



Woe " m c m 

Since [Pk,Pi] = {ih~e/c)(Ai,k — Ak,i) — (ihe/c)sikjHj, where Eikj is absolutely antisymmetric tensor, the 
property is proved. 

Lemma 3 The following relation holds: 

£o(s_2) 
2(m c) 



[HiM- = - J \ HxV)+p2(X,ExV) + 



^(p 3 divH + p 2 diyE + ip 3 {E,rotH)+ip 3 {E,TotE)) }. (B.6) 



[Wi.So]- = ' e °i g 2) — ^ p 3 [<S,H)<S,P)]_ + P2<S,S)<S,P>]- }>, 



Proof. Actually, 

r-H, . ,S*nl = - 

2too c m o c 

[<£,#>, <E,£)]_ = [(H,P)]-+i(Z,HxP)-i(Z,PxH) 
= 2i(Z,H x P) +ihdWH- h{Z,TotH), 

as was to be proved. 
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Lemma 4 The following relation holds: 



I^i'S] = ",1 J \ ( m oc)(-S xH- E Pl ) + P3 (E, E)P - P2 (S, ff)P - 



e (g - 2 ) 
2(m c) 2 

^grad(p 3 <£,I?>+/>2 <£,#>)!>. (B.7) 



Proof. Actually. 

2too c 



2m c 

Hi,—pifi 

m c 



ie (g - 2) 
2m c 



ip 2 [(E,H),P} + -ip 3 [(£,E),P] + } = 



' 1 1 '' 2) ip 3 (Z,E)fi-p 2 (Z,H)tj- l 4 grad(p 3 (£, E) - p 2 (t,H)) 



(to c) 2 |^-'-'- - — ' 2 

as was to be proved. 
Lemma 5 The following relation holds: 



nVw{§4[« ,s„]-}n +W = ^{(.-,E) + 7l f|> r ^/J)}. (B.8) 



Proof follows directly from (A. 9). 
Lemma 6 The following relation holds: 



+ 7 -^Hxff+—^— T (aJ}Hxf3\. (B.9) 
1 +7 J 

Proof follows directly from (A. 12) and (A. 9). 
Lemma 7 The following relation holds: 

n* + (t)[W!,5o]_n + (i) = - e A^^(aJ) r ^0,Hxf)- 

- ((J,H x P)^j - (a, (3 x (E x "P))| + 0(h). (B.10) 

Proof follows directly from (A. 7) and (A. 8). 
Lemma 8 The following relation holds: 

n* + (t)[Wi,^-n + (t) = ~ f {-m c((a, - 7-^ x 3 - 

-Y^-T^^^x^ - ^a,^)^ML-(/3,^)-(a,/3x#)^J+0(/l). (B.ll) 

Proof follows directly from (B.7), (A.7) and (A.8). 
Lemma 9 The following relations hold: 

nV(t)^ n+(t) = J-jV 1 ^) + (B.12) 

nVW^B+W =ff-{ff,fi—P + —P(aJ). (B.13) 

1 + 7 1 moc 

Proof follows directly from (A.9), (A.7) and (A.12). 



25 



Lemma 10 The vector 



satisfies the equation 



"' {(a, 0)E + Hxa) + ^ 2)l ((a, E) + (a, 0) 0, E) + 0,ax H)). (B.14) 



2t7JoC7 2tooc 
Proof. We express the vector ( in terms of a and obtain 

C = 3-—^- T (aJ) (B.15) 
1 + 7 

(see (A.l) and (A. 2)). Averaging (.9) with respect to and setting ft — > 0, we substitute (.15) into the 
relation obtained. Then 

= —^—Ua,0)E + Hxa}. 

m cj 



m c 



Similarly, we transform (.11) 

- m c((a, 0)E + H x a)^ 1 - 

- 7»W (7- 1 (a, P) + (a, E) + (a, 0) - (a, x H)^ = 

= e ^~ 2 \ l-\{a, 0}E + Hxa)+ 7 /3«a, E) + (a, (3) 0, E) - (8,0 x H))}. 
Substituting the latter relation into (2.30), we get (.14), as was to be proved. 
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